Observing the evolution of the CMB 
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Most cosmological parameters are expected to change significantly only on cosmological time 
scales, but given the large amount of information contained within the Cosmic Microwave Back- 
ground (CMB) sky, we can expect that changes in the CMB should be observable on much shorter 
time scales. Here we quantify this expectation, examining the detectability of the dominant effects 
on short time scales. We find that an ideal future experiment with currently achievable sensitivity 
could detect the changing dipole due to our galactic motion in about 10 years, but that it would 
take around 4000 years to detect a change in the higher order multipoles. 
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Introduction. One of the most scientifically profitable 
endeavors of recent times has been the increasingly ac- 
curate measurement of temperature anisotropics in the 
Cosmic Microwave Background (CMB) (see [lj, and ref- 
erences within, for reviews). These measurements have 
been at the cornerstone of constructing a standard model 
of cosmology, coined the ACDM model, which is seeded 
by near scale invariant adiabatic initial conditions cre- 
ated during inflation. 

The CMB sky is normally regarded as a "snapshot" of 
conditions on our last scattering surface (LSS) (together 
with some line-of-sight effects), with a finite information 
content which could be extracted by a sufficiently sen- 
sitive experiment. However, in recent work (2j, we have 
discussed extensively the time evolution of the CMB in 
the context of a ACDM cosmology. This raises the pos- 
sibility that one might learn more about the Universe by 
observing the changing CMB sky, i.e. by using the time 
domain, as is common with other astrophysical objects. 

The main evolutionary properties of the CMB can be 
classified into five categories: (1) an overall cooling of 
the mean temperature; (2) a changing dipole due to our 
evolving velocity with respect to the CMB frame; (3) new 
structure from probing a growing LSS; (4) geometrical 
effects of projection as the distance to the LSS grows; 
and (5) effects due to evolving structure along the line of 
sight. For small time scales (< 10 6 years), the second of 
these effects dominates on large angular scales, and the 
third on small scales !• As a result, if one were to take 
the difference between the observed CMB sky now and in 
the near future, the resulting map would be dominated 
by small scale features around the anisotropy damping 
scale of ~ 10 arcminutcs. 

These effects are to be distinguished from the be- 
haviour of the CMB at very late times in a ACDM cos- 



mology, which contains a future event horizon. Some- 
what akin to the stagnation a distant observer sees as 
objects approach the event horizon of a black hole, in 
the cosmological case an observer will see the pattern of 
CMB anisotropics "freeze in" (accompanied by ever in- 
creasing redshift) as the sphere of the LSS approaches a 
maximum radius on the de Sitter time scale (see, e.g., 

Since the sensitivity of CMB experiments has improved 
dramatically over the years, in this letter we extend our 
treatment in Ref. Q, and ask how long it would take 
to detect a change with current and future technology. 
This issue has also been briefly considered in Ref. Q , but 
here we take a more thorough approach to the detection 
statistics and also consider evolution of the CMB dipole. 

Detection criteria. The goal of this letter is to quan- 
tify when one could detect the temporal evolution of 
the CMB, so here we must first clarify this condition. 
A natural basis to expand the temperature fluctuations 
of the CMB are the spherical harmonics Yg m (fi), such 
that 5T(h, rj)/T(r]) = J2e m ae m (v) Y e m {h) , where ?] is the 
conformal time. The statistical properties governing the 
time evolution of the ai m coefficients are then encoded 
in the matrix (assuming isotropy) 
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where Cg V = Re (ai m {ri)a* em {ri')) s is a correlation func- 
tion. Here, (. . .)s denotes the ensemble average over ag m 
(signal) realizations. In particular, when ?/ = rj' this 
correlation function reduces to the ordinary anisotropy 
power spectrum, Cg V = C^rj). 

Wc find it convenient to define the power spectrum of 
the difference between a£ m modes at two times, Dg = 
(\Sa tm \ 2 ) Sl where 8a tm = a em (if) - a em (r]). Expanding 
this quantity, one finds De = Ce{r])+C£(r]')—2C^ V . Since 
Di contains covariance information between the a£ m s at 
different times, it should be better to measure this, rather 
than simply the change in power spectra, 8Ci = Ce(r]') — 
Ci(rj), if one were to attempt to detect a change. 
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Unfortunately, when making CMB sky maps to com- 
pare at different times, each of the a^ m modes will be 
corrupted by noise, such that a^ p = ag m + a^°l sc . Fur- 
thermore, since we only sample a finite number of ag m 
modes from the ensemble, one can only provide esti- 
mators of the true variance. If one constructs C" lap = 
l a £m P | 2 /(2^ + l)i then an appropriate estimate of the 
variance is C| st = C™ ap — Ne, where we assume the noise 
variance is also isotropic, i.e. Ne = (|a"„ se | 2 )Af- Here 
(. . .)jv refers to the average over the ensemble of noise re- 
alizations only. At this point, however, we need to make 
the distinction between two situations involving C| st : 

(A) The averaging is taken over all possible signal and 
noise realizations, and is denoted (...). The ensemble 
average of the variance estimate is (C£ st ) = Cg, and so 
the estimate is unbiased. In standard parameter estima- 
tion analysis, one is often interested in the width of this 
distribution, i.e. (ACe) 2 = ((C c e st - Ce) 2 ). In a full sky 
CMB experiment it can be shown that this quantity is 

given by [!, 0] ACe = J jeTi^e + Ne) . Hence, even if 

Ni is zero there is still sampling uncertainty of the un- 
derlying signal variance. This is commonly referred to as 
"cosmic variance" . 

(B) In the second situation the signal component is 
fixed and we average only over the ensemble of noise re- 
alizations. Hence, we are interested solely in the uncer- 
tainty of the measurement of a given signal due to noise 
corruption, not the uncertainty in its underlying value. 
To emphasize the fact that we now regard the signal as 
fixed, we denote the realization of the ai m component 
that we are dealt by ai m . If we construct the variance 
of this component by Ce = ^2 m \ae m \ 2 /(2£ + 1), then we 
now find that (C£ st )jv = Cg. Furthermore, the width of 
this distribution, (ACf ) 2 = ((Cf * - C e ) 2 ) N , is given by 

ACf = ] J™!L(2C t + N t ). (2) 

In this case, the uncertainty in measurement vanishes 
when the experimental noise is zero, as one would expect. 

Applying these situations to the time evolution of the 
CMB, one can therefore ask two different questions: has 
the realization of the sky that we have been dealt changed 
within experimental noise uncertainty (corresponding to 
case B), or has our realization changed to such an extent 
that we can conclude that the ensemble averaged spec- 
trum has changed (corresponding to case A)? Clearly, the 
former criterion is easier to meet than the latter. 

To illustrate the size of the errors given by ACe and 
AC^ , in Fig. Q] we plot these quantities along with 
the primary CMB temperature anisotropics for an in- 
strument with Ni = 1 x 10~ 4 fiK 2 /T 2 . Such a noise 
level will be obtainable by the upcoming Planck exper- 
iment 0- On large scales, ACe is comparable to the 
primary anisotropics due to sample variance. However, 
even though we may be uncertain about the underlying 
variance, we can in principle measure Ce with high pre- 
cision: AC^ is considerably lower than ACe. 
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FIG. 1: Primary CMB temperature anisotropies (solid curve) 
compared to the error in Ce- The dashed curve shows ACe, 
and the dotted curve AC? . We assume Planck noise levels. 



Equivalent expressions apply in pixel space — in this 
case the temperature at each pixel i of the CMB sky map 
contains a signal and noise contribution: (6Ti/T) mlip = 
(STi/T) + (STJT) 1 ™ 30 . The expected signal variance 
over the ensemble for each pixel is Cs = ((STi/T) 2 ) g = 
J2e (2^ + l)CV/(47r) (i.e. the total power in the map), and 
the noise variance C N = ((5T™ iec /T) 2 ) N = cr 2 ix . If one 
constructs C map = Ei^T™* / T ) 2 /N pix , where N pix is 
the total number of pixels in the map, then an appropri- 
ate estimate of the variance is C s st = C" nap — Cm- If one 
defines Cs = E^S^/T) 2 /N p [ x , then averaging over the 
ensemble of noise realizations gives (C§ st )jv = Cg. The 
equivalent error on the signal variance to that given above 

becomes ACg = ^2C N (2C~g + C N ) /N pix . Analogous 
quantities arc defined for the difference map power D. 

Using these expressions we can estimate the time-scale 
for detectability. However, it will first be instructive to 
consider a special case before considering higher order 
CMB modes. 

Time varying dipole. The CMB dipolc (£ = 1) is deter- 
mined by the Doppler shift of the approximately isotropic 
CMB radiation field due to our observer motion with re- 
spect to its rest frame. The magnitude of the dipole 
variance is much higher than modes with I > 2, and the 
evolution has contributions on much shorter time scales. 
Indeed, the time variation of the Earth around the Sun 
has already been well measured and is used to cali- 
brate CMB experiments, so here we consider larger scale 
solar motions. 

For simplicity, we decompose our observer velocity into 
two components: v = v Slm -GC + vqc-cmb- The first 
component is due to our solar motion around the Galac- 
tic centre (GC), and the second is due to motion of the 
GC with respect to the local group and other large scale 
structures. We assume our galactic motion is confined to 
the plane, with a tangential velocity of 222 kms -1 at a 
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radius of 8.5 kpc, giving a galactic rotation period T g of 
2.35 x 10 8 years. We assume the v GC _ CMB component 
remains constant on galactic time-scales. 

The best estimate of the measured a\ m components 
come from the WMAP satellite @, and are given 
in fractional temperature units by (ai-i, aio, an) = 
(v y ,v z ,v x ) = (-816.0,919.3,-87.8) x 10~ 6 , with 
around 0.5% calibration uncertainty. The predicted evo- 
lution of the dipole components due to galactic motion 
is then (oi_i, aio, a n) = 740.5 (cosx, 0, sinx) x 10~ 6 + 
(-1556.5, 919.3, -87.8) x lO" 6 , where x = 2i:t/T g and 
t is the local observer time. 

For small time periods t <C T g , we can compute the 
evolution in the dipole variance C\ = J2 m a i™| 2 /3 to be 
d{t) w Ci(t ) - 1.16 x 10- 15 £/years, where Ci(i ) = 
5.062 x 10 -7 . The evolution in the difference power is 
Di = 1.30 x 10~ 22 (t/years) 2 . Although SCi is larger 
than D\ over small time-scales, the larger noise error 
means it is a less sensitive statistic, as we show below. 

The simplest estimate of the time-scale for detectabil- 
ity is as follows. We assume that we have a full sky data- 
set with no noise today, and use the predicted evolution 
of the ai m coefficients to generate another noiseless sky 
at some time in the future. Assuming the pixel noise 
is uncorrelated, one constructs the hypothesis that the 
two skies are different within the expected noise. So we 

form x 2 = Ef pix [r i /T(t)-T l /T(t )] 2 /(2C 7V ) ! where 
the factor 2 in the noise arises since wc arc differencing 
maps. Under the null hypothesis, the expectation for this 
quantity is zero and the standard deviation is y / 2N p - lx , 
so one can evaluate the significance level a that the two 
skies are different. Note that for the dipole, we only have 
noise averaging, as the signal is deterministic, unlike the 
£ > 2 modes. 

The noise contributes on all angular scales (not just 
the £ = 1 mode), and may swamp the signal variance of 
the difference map, D$- In the noise dominated regime, 
the simple % 2 statistic is equivalent to the signal-to-noise 
(S/N) ratio a = Ds/ADg. In an experiment such as 
Planck, the large number of pixels («4x 10 6 ) will beat 
down the noise level, and one finds a rj 2 after t = 2 x 10 
years, even though the map is noise dominated. 

Much of this noise, however, resides at £ > 2. If one 
performs a harmonic transform on the noisy map, one 
can recover a higher S/N by isolating the £ = 1 mode. 
In this case, the uncertainty in the signal variance given 
the noise is AD% = 3ADf /(4tt). At t = 2 x 10 4 years, 
by filtering small scale noise, the detection significance 
increases to a w 40. 

Finally, one could also compare the change in power 
SCi. The uncertainty in Ci due to noise at each time 
is ACf . Since Ci(f) « Ci{t ), ACf (t) ra ACf (to) = 
ACf, and the S/N is a = |<5C*i |/(2ACf). After t = 
2 x 10 4 years, an experiment with Planck noise gives 
a ~ 4. Therefore, the loss of time covariance information 
has led to a reduction in statistical significance over the 
filtered difference map method. If one instead assumed 
the signal was stochastic, and used ACi instead of ACf , 
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FIG. 2: Ideal case detection signal-to-noise for observing a 
change in the CMB dipole. Thin lines show an instrument 
with Planck sensitivity, and thick lines an instrument 100 
times as sensitive. Solid curves show the significance using 
the noise dominated difference map, and the dashed curves 
after removing small scale noise. Dotted curves are for the 
dipole magnitude alone, and hence contain no covariance in- 
formation. 



one would obtain a tiny significance, as one is estimat- 
ing the probability that two uncorrelated samples have a 
different underlying ensemble variance. 

In this idealized case, we plot the detection S/N in 
Fig. [2] for Planck noise levels, along with an instrument 
100 times as sensitive (this is current generation technol- 
ogy [l(|). In principle, such an instrument could observe 
the evolving galactic dipole on a ~ 10 year time-scale. 

An important systematic in any such experiment 
would be calibration uncertainty, the accuracy required 
being many orders of magnitude better than what has 
been achieved so far. There are two contributions to 
dipole evolution — a change in amplitude due to changing 
speed relative to the CMB frame, and a change in orien- 
tation due to rotation. By normalizing Ci(t) = Ci(to), 
hopefully eliminating calibration issues, the only con- 
tribution is due to rotation. Due to our position in 
the galaxy (our velocity is mainly against the bulk ve- 
locity flow), the rotation component actually dominates 
currently, and amplitude normalization only leads to a 
degradation in S/N of < 0.5%. Another more insidi- 
ous effect is time-varying foreground emission. As well 
as the motion of Galactic clouds, the observed fore- 
ground pattern will change on the same time-scale as 
the CMB dipole. It is difficult to place a degradation 
factor on these effects, although one would hope that 
multi-wavelength observations would minimize them. 

Higher order anisotropics. In contrast with the dipole, 
the evolution of modes with £ > 2 depends only on the 
set of cosmological parameters. For small time intervals 
t (< 10 6 years), the spectral shape of the change in vari- 
ance SCi and difference power Dg are fixed, with the 
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t (years) 

FIG. 3: Ideal case detection signal-to-noise for observing a 
change in the higher order [l > 2) modes. Labels are the 
same as in Fig. [21 

amplitudes scaling oc t and oc t 2 respectively The 
shape of Dg is relatively featureless, being approximately 
constant up to the anisotropy damping scale. On scales 
smaller than this (I > 1000) the power falls off rapidly. 

Truncating the spectrum at I = 2000 (the approxi- 
mate resolution of Planck), the total expected difference 
power is given by D s = 3.11 x 10~ 23 (i/19.8years) 2 . The 
prefactor in this expression is the dipole signal variance, 
3-Di/(47r), so it would take about 20 times longer than for 
the dipole to obtain the same S/N a = Ds/ADg. This 
is borne out in Fig. [3J with around 4000 years required 
to detect a change in the I > 2 modes at the a = 2 level 
using current detector technology. 

In the case of the dipole, we were able to reduce this 
nominal value significantly by removing small scale noise. 
We follow a similar procedure here, computing the ex- 
pected S/N at each £ mode, ag = Dg/AD^ , and dis- 
carding modes with low significance. In order to com- 
pute the combined significance, we first rank the modes 
in terms of ai. Applying a weight to the signal at each 
I mode, given by wi = (ADf)~ 2 , the combined signal 
and noise components are then given by wiDi / £\ wi 
and (52iW()~ l t 2 respectively. By iterating through the 
modes, starting from high at, we compute the com- 
bined significance until the maximal S/N is obtained. 
Typically, this procedure retains modes in the range 
20 < I < 500, so the limiting factor in any such mea- 



surement is instrument sensitivity rather than resolution. 
Unfortunately, since the signal contributes on a range of 
scales, similar improvements to the dipole are not possi- 
ble here, and the increase in S/N is small, as shown in 
Fig. [3j Our results for I > 2 are consistent with Ref. [H , 
who use ADi as the noise, without the cosmic variance 
term. Since the signal at each i mode is noise dominated, 
ADg m ADf in this regime. 

Following the same method using the expected change 
in variance, 5Ce, such that the S/N at each i mode is 
oil = 5Cg /(2AC e N ), the time for detection increases by 
around a factor of ten before one could detect a statistical 
change. 

Where does the stochastic nature of the signal for the 
£ > 2 modes enter into this treatment? For small time 
intervals, the evolution of anisotropics is primarily driven 
by probing new structure at the surface of last scattering. 
The spectrum Di of this structure has a corresponding 
sample variance, equal to 2D 2 / (21+1) at each t mode. In 
these estimates, we have assumed that the future signal 
we would measure is at the peak of the likelihood function 
(equal to Di). The uncertainty in the exact signal we 
would measure induces a corresponding spread on the 
S/N, but we found this to be negligible after performing 
realizations of the future sky. 

Conclusions. The time domain provides critical infor- 
mation in many branches of astrophysics, even cosmology 
(distant supcrnovae and active galactic nuclei being ob- 
vious examples). However, the Universe as a whole is 
evolving slowly and so one expects this evolution to be 
difficult to observe, as is the case with the Sandage-Loeb 
test [ill [HI for the increase in redshift of fixed sources. 

In this letter we answered the question of how long one 
would have to wait to see new information from the gold- 
mine of cosmological data that is the CMB. The large 
quantity of data available means that we could observe 
a change on time scales much shorter than cosmological. 
By considering optimal estimates for differences in ob- 
served skies we estimate that, using currently available 
detector technology, the dipole can be observed to change 
due to our Galactic motion in about a decade, whereas it 
would take thousands of years to detect even a statistical 
change in the higher order multipoles. 
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